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' Abstract 

^ ; 

Q-i| We present the results of the quantum field theory approach to nuclear Yukawa model obtained 

■ by standard lattice techniques. We have considered the simplest case of two identical fermions 
interacting via a scalar meson exchange. Calculations have been performed using Wilson fermions 
in the quenched approximation. We found the existence of a critical coupling constant above which 
the model cannot be numerically solved. The range of the accessible coupling constants is below 
the threshold value for producing two-body bound states. Two-body scattering lengths have been 

^ I obtained and compared to the non relativistic results. 

^ ! 1 Introduction 

The application of lattice techniques to nuclear physics is nowadays an active and fruitful field of 
^ \ investigation. This activity has been extensively developed in the last ten years and covers different 

■ aspects of the problem. 
The first one aims to obtain relevant nuclear properties from ab-initio lattice QCD calculations 

(LQCD). In this approach, the elementary fields are quarks and gluons and all the numerical simu- 
^ ' lations depend on the very few QCD parameters: the bare coupling constant and the quark masses. 

■ The first task there is to generate the nucleon from its elementary constituents and thus this approach 
is necessarily limited to very light nuclei. First unquenched results on NN scattering length [T] and 
NN potentials [21 [3] have been obtained although still for large pion masses. Very recently, the bind- 
ing energy of the A=3 and A=4 nuclei have been computed [1]. Although being performed in the 

^ . quenched approximation, this is a result that seemed out of range just few years ago. 
^ \ A second approach has been developed using the same techniques but in the framework of effective 

field theories, i.e. using nucleon and mesons as elementary fields. Lattice effective field theory (LEFT) 
has been first applied to study nuclear [5] and neutron matter [6] and was latter adapted to light 
nuclei [7]. The action describes non relativistic nucleons interacting via realistic, chiral inspired, NN 
potentials. These potentials depend on a large number of parameters depending on the order of chiral 
FT and have proved to be very successful in the non relativistic Faddeev-Yakubovski description of 
the A=2,3,4 nucleon problem. LEFT is a quantum mechanical description of a many body system, 
equivalent to a ladder potential models, that take the simplicity of the lattice techniques to go beyond 
few-body methods. At present this allows to treat bound states of nuclei up to A=12 nucleons and 
can be extended well beyond. A review on this approach can be found in [8]. 

Our aim in this work was to consider the simplest Quantum Field Theory (QFT) model of inter- 
acting fermions which was at the origin of nuclear forces: the Yukawa model. Since Yukawa pioneer 
work [9], the meson-exchange interactions constitute the starting point for building the NN poten- 
tials [ini [HI [12] which, inserted in Schrodinger-like equations, provides an "ab-initio" description of 
light nuclei up to A ~ 10 [13]. The potential approach, however, takes into account only a small, 
though infinite, fraction of diagrams of the perturbative series ~ the ladder sum. This represents 
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a severe restriction of the interaction, specially taking into account the large values of the coupling 
constants involved. Chiral inspired NN models |14[ I15[ [IB], which can be formally distinguished from 
the traditional meson-exchange ones, suffer from the same restrictions. 

To incorporate the full content of the meson-exchange Lagrangian we have used the standard lattice 
techniques [T7], developed in the context of QCD. They are based on a discrete Feynman path integral 
formulation of QFT and provide nowadays a genuine way to solve non perturbatively such problems. 
Preliminary results of this work can be found in \18\ [2T\ 122] . A similar study was undertaken in [23] 
in the frame of a purely scalar c/i'^x model. 

A few lattice investigations of the Yukawa model with an additional Xcp'^ term have also been 
performed some time ago in the spontaneously broken phase |24j . Similar models were studied in the 
past (see for instance [T7] and references therein), and the existence of numerical instabilities beyond 
a critical value of the coupling were found. Different schemes were used for the discretization of the 
meson field, obtaining critical values of the lattice coupling slightly different, but corresponding to the 
same continuum limit. More recently the phase structure of a chirally invariant lattice Higgs- Yukawa 
model was studied to establish Higgs boson mass bounds [251 ES]. None of these studies was however 
interested to investigate the existence of two-fermion bound states in the original Yukawa model. 

The plan of the paper is as follows. In section [2] we describe the model and its discretization as 
well as the approximations in use. Section [3] contains the detail of the Monte Carlo simulations and 
the results for the renormalized fermion mass. The two-body binding energies and scattering lengths 
are presented in section U] and compared to the non relativistic results. Some concluding remarks are 
finally drawn in section [5l 

2 The model 

We have considered the simplest renormalizable Quantum Field Theory describing a fermion ('J/) - 
meson (^>) interaction. It is given by the Lagrangian density: 

C = Cd{^, ^) + Ckg{^) + Cii^, ^, $) , (1) 

with, in Euclidean space, 

Cd{^,^) = ^(a^T'^ + mo)* , (2) 

^KG^ = ^ + , (3) 

Ci{^, ^, «>) = c/o"^r$^ + Ao$^ , (4) 

The fermion field is supposed to describe a nucleon (N) and the meson field a - more or less fictitious 
- scalar particle (a) responsible for the attractive part of the NN potentials. The Lagrangian depends 
on four parameters: the fermion niQ and meson fiQ masses and two dimensionless coupling constants 
go and Aq. The Yukawa coupling in Eq. @ admits several possibilities depending on the choice of 
the r matrix: P = / for the scalar case and P = ij^ for the pseudo-scalar one. In this work we will 
restrict to the scalar coupling, for it is known to produce a stronger fermion-fermion attraction. 

This theory can be treated perturbatively by computing order by order the contributions in qq, as 
is done for example in QED for computing the anomalous magnetic moment of the electron. Bound 
states nevertheless appear only non-perturbatively, when contributions at all orders in the coupling 
are taken into account. 

Non-perturbative tools are extensively used in the context of QCD. Among them, the lattice 
techniques provide a reliable numerical method for solving any QFT on a discretized Euclidean space- 
time. A very basic description of these techniques is given below. The interested reader can found a 
detailed explanation in e.g. [T7] . 
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The vacuum expectation values of the operators involved are obtained in the Path Integral approach 
which consist in computing integrals like: 



(0(^,^,$)> = ^ y"[d^'][(i'^][(i$]0(^','f,$) e-^^[*'*'*] , (5) 

where the Euclidean action plays the role of a probability distribution in a Monte Carlo simulation. 
The discretized Euclidean action Se can be written, according to Eq. ([1]), in the form: 

Se = ^ C = Sd + Skg + Si 

X 

where a is the lattice spacing and x denotes a point with coordinates = an^ {^l = 1,2,3,4) and 
riij, = 1 . . . , L^. In practice we have taken equal spatial dimensions Li = L2 = = L and a temporal 
one with L4 = 2L. 

For the free Dirac action Sd, we have used Wilson fermions. They consist in adding to the, naively 
discretized, derivative term of the Dirac Lagrangian ([2]) a Laplacian operator in order to remove the 
spurious poles at the boundaries of the Brillouin zone [T7] . The free Dirac action is then written as a 

I 3 

bilinear form in the dimensionless fermion fields ip = \ 

SD=^i^xDYy^y (6) 

xy 

where 

DZ = ^x,y - K ^ [(1 - 7^) (5^,j;„^ + (1 + 7m) ^x,y+tj\ (7) 
A* 

is the Dirac- Wilson operator, and 

K = ^ (8) 

8 + 2amo 

is the hopping parameter. 

In terms of the dimensionless meson field 6 = a$, the discrete Klein-Gordon action reads: 



(9) 
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X 

The interaction term takes the form 

Si = Tpx(l)xil^x + Ao ^ (/>^ (10) 

X X 

where = '^ngo is the lattice Yukawa coupling. 

Taking together the fermion meson ([9|) and interaction (jlOp terms, the lattice Euclidean action 
can be written in the form: 

Ssilp^lpA) = '^'4'xDxy'tpy + Sm{(P) (11) 

xy 

where Sm includes both the Klein-Gordon and the X(j)'^ term 

SM = SKG + XoY.^t (12) 
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and the fermionic part is written in terms of the interacting Dirac operator 



(13) 



Notice that the model is now made dimensionless. When needed, the physical quantities ~ masses, 
energies, etc. - will be given in terms of the lattice spacing. 

One of the most demanding issues when computing the vacuum expectation values ([5]) comes from 
the Grassmannian character of fermion fields, that have to be integrated out by algebraic methods. 
For example the fermion propagator, corresponding to 0{'tlj,'ip,(p) = tpxi^y, 



implies the evaluation of a determinant and inverse of an operator that, even for moderate lattices, 
V ~ 24^, has a dimension of ~ 10^. Moreover, if a Monte Carlo simulation is to be done using Eq. 
(|14p . the probability distribution for meson configurations is given by e~'^A/(<A)-iog(det(_D))^ what means 
evaluating a large determinant in every Monte Carlo step . This can be avoided by the use of Hybrid 
Monte Carlo techniques that nevertheless are the main source of time spent in the simulation. This 
task is considerably simplified in the "quenched" approximation that, from the computational point 
of view consists in setting det(D) independent of the meson field in the fermionic integral. 

From a physical point of view, the quenched approximation avoids the possibility for a meson to 
create a virtual nucleon-antinucleon pair (p — >• ipip (see figure [1]) . Due to the heaviness of the nucleon 
with respect to the exchanged meson this approximation is fully justified in low energy nuclear physics 
and implicitly assumed in all the potential models. 



Figure 1: The quenched approximation neglects the possibility for a meson <I> to create a virtual 
fermion-antifermion pair 

We have furthermore chosen in our simulation to neglect the meson self-interaction term setting 
Ao = in (|12p . This choice is consistent only in the quenched approximation. In a full QFT treatment 
of this model, the fermions loops will generate meson self- interactions, that require a Ao</>'^ counter-term 
for renormalizability. 

The model depends on three dimensionless parameters go, afio, and aniQ usually set via the hopping 
parameter These parameters appearing in the Lagrangian are not physical: they are modified by 
the interaction so that they have to be renormalized. Our first task to map the bare quantities into 
the renormalized parameter space 



In the quenched approximation and for Aq = 0, mesons do not interact each other and therefore 
its mass renormalize trivially a/x/j = a/io- In the following we will omit subscripts and denote the 
dimensionless meson mass by a/x. The lattice spacing a does never appear explicitly and it has to 
be fixed by setting a physical dimensional quantity. We do that by identifying the unchanged meson 
mass /i to a physical meson of /i = 0.65 GeV, a typical value used in the NN models. If we are using in 
our simulations a/i = 0.1, the lattice spacing is given by a = ^ = q gs gcV ^ ^'^^ GeV~^ ~ 0.03 fm. 

Coupling constant renormalization is a more involved issue. Renormalized coupling constants were 
computed in a previous work [18]. For the scalar coupling no sizable effect of renomalization was found. 




(14) 




(a/io, 50, arno) =^ {afiR, qr, amn) 
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i.e. gn ^ go in a wide range of momenta. From now on, we will denote by g this unique coupling 
constant. 

The remaining task for finding the adequate parameter space is thus to compute the renormalized 
fermion mass as a function of the bare parameters. This will be described in the next section. 

3 Monte-Carlo simulation 



3.1 Generating field configuration 

In the above defined conditions, the generation of meson field configurations is straightforward. This 
can be seen by writing the Klein-Gordon action ([9]) in Fourier space 



_ 2 



Skg = 



(15) 



where d>k denotes the Fourier transform of the meson field 



-ik-x 



and we have introduced the quantity 



1 



which depends on the lattice discretized momenta 

k,, , 27r 



k^ = 2 sin '-^ k^ 



One can see from equation (jl5p that the different Fourier modes (pk of the meson field appearing in 
the action are decoupled and can be generated independently. 

The Monte Carlo algorithm becomes then trivial as it is enough to produce at each lattice point 
k, independent complex scalar fields with a probability density given by 



exp 



1 



2 



exp 



lRe[^ 
2 



exp 



l lm[0i] 
2 c^l 



(16) 



i.e. centered Gaussian distributions, both for their real and imaginary parts, with a variance (7^ 
depending on k and the constraint imposed by the reality of <j)x ■ This method generates configurations 
that are statistically independent, thus saving a large amount of computing time with respect to the 
Metropolis algorithm. 

The scalar fields in configuration space (f)x are finally obtained by performing an inverse Fourier 
transform on 4>k- It follows from the particular form (I16p that (px are also centered Gaussian with a 
width a, independent of x, given by 
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V ^ k'^ + a? 



(17) 



Note however that the (j)x are now correlated. The correlation function is given by 

1 

^ k"^ + a^/i^ 



A(x - y) =< 4>x4>y >= Y 



Jk(x-y) 



(18) 
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which is the scalar propagator in configuration space. It is interesting to note that the parameter a 
in Eq. (|17|) is related to the lattice regularization of the potential at the origin, 

y(0) = VA(0) = -gW . (19) 

The a/i-dependence of a is displayed in figure [2] for different values of the lattice size, L. It behaves 
like 1 /afi in the two trivial limits a// — )■ and a// ^ 1 with a plateau in between, which is the region 
we are interested in. The parameter a allows a discussion of the discretization and finite volume errors 
in terms of physically well defined limiting cases: 

• For large values of a/i, the k"^ term in (I17p becomes negligible and <7 = ^ for any value of L. 
This correspond to a contact interaction between fermions. 



• In the limit a/x — t- the sum p7|) is dominated by the mode 4'k=o which generates a behavior 




Figure 3: Discrete Yukawa potential (VLatUce) 
Figure 2: a/i-dependence of the variance for the for L = 48 and a/x = 0.1 (Black dots) compared 
meson fields in configuration space. to the continuum one (Solid line). 

The non-relativistic Yukawa potential can be computed for the discrete lattice using Eq. (llSp but 
summing only over spatial directions. The resulting potential {ViatUce) for afi = 0.1 and L = 48 is 
represented in figure [3] (Black dots) and compared to the continuum result (Solid line). The main 
effect of discretization is the regularization at the origin and is seen to be negligible beyond the very 
first points. This potential will be used in section |4] to compare continuum and lattice results. 

The appearance of volume effects depend crucially on the afi value: they are very small for large 
values of afi but important when afj. — )• 0. As a matter of fact for a given value of a/j,, there is a 
minimal lattice size L below which the lattice artifacts are dominant. We took this constraint into 
account in the present work. 

3.2 Zero modes of Wilson-Dirac Operator 

When computing physical observables, the integration over the fermionic fields is performed analyti- 
cally and the result is expressed in terms of the inverse Dirac operator (jl3p . This has been explicitly 
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done in eq. ()14p for the fermion propagator which constitutes the building block of the lattice sim- 
ulations. When working in the quenched approximation, one can set det(D) = 1 and the relevant 
numerical task is thus reduced to compute D~y [0] for an statistical ensemble of meson field configu- 
rations. 

1 r 1 ^ 

S{x,y) = ^ m D-yH4>) e-'^'^'^^ « -Y^D-y\^i) (20) 

1=1 

Due to translational invariance one is left in practice to compute 5[(/)](x,0) = D~q^ [(/>], that is to 
solve the linear system: 

D,M)Sx{^) = 5,0 (21) 

It is worth noticing that in the full QFT formulation every configuration is weighted by the 
determinant of the Dirac operator D and therefore the configurations yielding an ill-conditioned linear 
system ()2ip . i.e with det(D) ~ 0, do not contribute to the functional integral. In the quenched 
approximation, however, this is no longer true and "ill-conditioned configurations" can be sampled. 

As a practical measure of the "ill-conditioness" of D we have considered its "condition number" 
defined as the ratio between the largest to the lowest eigenvalue modulus [E]. The largest is this 
number the more difficult is to solve the linear system. Depending on the method used for that 
purpose, either the algorithm cannot find the solution, or the round-off errors make the solution 
wrong. In exact arithmetic the condition number measures how the solution changes when the second 
member of linear system slightly changes. 

We have found that such "ill-conditioned configurations" appear in the Yukawa model for almost 
any k when gi > 0.6. In this case the inversion of the Dirac operator becomes in practice impossible 
|20j . For illustrative purposes, we have plotted in figure S] the condition number of D as a function 
of the lattice coupling constant for an ensemble of L = 8 configurations at fixed value of k. As 
one can see, the condition number of a given configuration diverges on a discrete set of gi values for 
9l ^ 0.6 indicating the practical impossibility to compute the nucleon propagator. The precise gi 
values where this divergence occurs depend on the particular configuration, on the values of k and afj, 
and on the lattice size. It turns out however that the situation described in figure U] is generic for the 
quenched Yukawa model. 
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Figure 4: Conditioning number as a function of gL for a fixed value of k = 0.11 and V = 8'^ and for 9 
different meson field configurations. 
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The existence of zero modes in the quenched Yukawa model was already found in [23 but the 
numerical results performed with very small lattice let these authors suspect the existence of a second 
allowed region at large values oi gi. It can be shown however that this result is a volume artifact as 
the second region disappears exponentially with the number of lattice sites in the thermodynamical 
limit. This issue as well as some properties of the interacting Dirac operator (jl3p . in particular the 
fact that its spectral properties do not depend independently of n and g but rather on their ratio, will 
be treated in detail in a forthcoming publication [20j . 

As a conclusion, the numerical simulations in the quenched Yukawa model are limited to values 
of the lattice coupling constant gL ^ 0.6. Using a typical value of k = 0.1, this corresponds to 

2 

g' = 1^ < 3, that is a = < 0.7 which is of the same order than the aqcD in the nonperturbative 
region. 

3.3 Renormalized Fermion mass 

Renormalized particle masses are obtained in Lattice QFT by considering the time evolution of the 
correlator matrix defined as: 

C{p,t) = X;(j(x)J^(0))e^?^- , (22) 

X 

where creates a particle state at the origin and J destroys it at x. The tensorial indexes - depending 
on the type of particle - are implicit and the vacuum expectation value (• • •) is obtained through an 
average over field configurations. It can be shown that the correlator matrix ()22p has contributions 
from all the particle states n satisfying (0| J|n) / 0, and has the form 



T:C(p,t) =^c„cosh£;„U-|j 



(23) 



For p = 0, it behaves as a sum of hyperbolic cosine with the rest mass rrii of the particle states: 



Tr C{t) = Co cosh arriQ ( * ~ j + "^i ^^sh. ami f ^ ~ ) + ' ' ' (^'^) 



Using the above equation at two consecutive times one can extract an effective mass 

C{t) 



anicsit) = T 



C(t + 1) 



(25) 



which, at large enough euclidean times, will display a plateau region that will be identified to amo 
value. 

In the case of one fermion state ( J = ^'), the correlator matrix is he 4 x 4 matrix: 

C{t) = Y,S{x,0) (26) 

X 

where S is the propagator defined in Eq. [TH The case of two fermion states will be discussed in 
section [H 

The fermion masses extracted in this way are presented in figure [5] for afj, = 0.1, a lattice size 
L = 24 and several values of the lattice coupling ranging from g^ = 0.0 to gi = 0.5. In the free case, 
this mass is already different from the bare one aniQ due to lattice artifacts and it is given by: 

am„ = log (l + 1(1-1)) (27) 
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with Z = 1 and Kc = 1/8. This expression is used to fit the interacting masses in terms of two 
parameters Z[g) and Hc{g)- The result of this parametrization is indicated by full lines in figure O 
Note that at k = Kc the renormalized fermion mass vanishes. As one can see, renormalized masses are 
smaller with growing values of the scalar coupling constant, or in other terms, Kc{g) < «^c(0). This 
indicates that the renormalized nucleon mass is made lighter by a scalar coupling. 

The fitted coefficient Ate (5) is presented in figure O This coefficient can be calculated in lattice 
perturbation theory which provides a test of numerical simulations. It is quadratic in the coupling 
constant ^ 

Kc{g, L, afi) = -- C2{L, a^j)g^ + o(/) (28) 

with coefficient C2 depending on the scalar mass o/i and on the lattice size L. 




Figure 5: Fermion mass vs the hopping param- 
eter for several values of Yukawa coupling aver- 
aged over 100 samples each one. The line is a 
best fit according to Eq. (j27|) . 



Figure 6: Values of Hc obtained with L = 24 
and a/x = 0.1 extracted from data in figure [5] 
according to Eq. ()27p . The line represents the 
perturbative behavior given by Eq. ([28]) . 



The region of parameter space to perform the numerical simulations of physical interest is limited 
by the constraint amji > which corresponds to k G [0, Kc{g, L,afi)]. In order to avoid large lattice 
artifacts the condition am^ <^ 1 must hold, what restricts the useful parameter space to a narrow 
band in the vicinity of Kc with the constraint indicated in the preceding section, say gL ^ 0.6. 



4 Two-fermion states 



Our main interest in this work is to study the interactions between fermions and thus the properties 
of the two-body states. In this respect, our reference will be the results provided by potential models 
in the non relativistic ladder approximations. These are summarized in what follows. 



4.1 Non relativistic results 

Let us first consider the non relativistic system of two particles with equal mass m, interacting by a 
Yukawa potential of strength g and range parameter /i 

V{r) = -^ — 
4tt r 
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The binding energy (B) and scattering length (ag) are given by 

2 



B 

ao 



m 



e(G) 



1 



KG) 



(29) 
(30) 



where e(G) and A(G) are respectively the binding energy and scattering length of the dimensionless 
S-wave Schrodinger equation. 



-e + G 



u{x) = 



(31) 



with a coupling constant G related to the original parameters {m,g,fj.) by 

G 



m 
47r /i 



The functions e(G) and \{G) are displayed in figures [7] and [HI The convention used for the scattering 
length corresponds to 5{k) = — oo/c + o{k'^). The critical value for the appearance of the ground state 
is Go ~ 1.680. At this value A(G) has a pole and it can be shown that for small values of G one has 



A(G) = -G + o(G2) 
which corresponds to the Born approximation. 



(32) 





Figure 7: Ground state binding energy of the di- 
mensionless non relativistic Yukawa model (j3ip 
as a function of the coupling constant G. The 
appearance of the first bound state corresponds 
to Go = 1.680. 



Figure 8: Scattering length in the dimensionless 
non relativistic Yukawa model as a function of 
the coupling constant G (solid line). The Born 
approximation is indicated by the dashed line. 
The singularity corresponds to the appearance 
of the first bound state. 



4.2 Binding energies 

We will restrict ourselves to study the system of two identical fermions, say A^A^, in the J'^ = O"*" state. 
The general form of the interpolating field J{x) for a two fermion state reads 

J{x) = ^ ilJa{x)TapilJii{x) (33) 

a/3 
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where T depends on the quantum number of the state. For a 0"*" state one has F = ^727075- In the 
matrix form l\3'6\i can be written as 

J{x) = i^iVi;^ 

Using the interpolating field (j33p . the NN correlator becomes 

C{t) = ^(0|J(0)jt(x)|0) = ^(0|VorVoVix7or7oV^x|0) = 2 Tr^(0|FS(x,0)FS(x,0)|0) 

XX X 

where S{x, 0) denotes the fermion propagator. 

In order to decrease the contribution of the excited states we have introduced a smearing procedure. 
This consists in modifying the interpolating field J{x) in the following way: 

J{x) = ^fiR)^,,tT^P^^n,t (34) 
R 

where f{R) is some smearing function to account for the spatial extension of the state. As we are 
interested in the s-wave, the smearing can be done over the whole timeslice in an efficient and inex- 
pensive way by choosing the smearing function to be a constant, f{R) = After performing 
the appropriate Wick contractions, the time-correlator results: 

C{t) = (TCi{t)TCi{t)-TCi{t)T'Cl{t)^ (35) 

where the 4x4 matrices Ci(t) and Ci{t) are respectively: 

at) = (36) 

w 

Ci{t) = J^5(f,t) (37) 

w 

and S{x, t) is the solution of the linear system 

D-SS^,^ = j^Y.^'''^R^oO- (38) 
R 

This three-dimensional smearing efficiently removes higher energy contributions to the correlator 
as it can be seen in figure [9] where the effective mass (I25p is plotted both for the local and smeared 
interpolating fields. It can be seen that for large t values the effective mass tends to a plateau which 
actually defines the mass of the state. If the local sources were used we would have needed larger 
values of t to be able to find a plateau in the masses. 

An interesting property of this smearing is that the two-body free correlator is the square of the 
one-body one at any value of t, i.e., the effective mass of the two-body state is rigorously constant and 
exactly twice that of one particle. 

The binding energy of a two-fermion state with total mass m2 is defined as i? = 2mi — m2 where 
nil denotes the mass of the fermion. In figure [10] we show this binding energy as a function of the 
lattice size La/j, for a given set of parameters. The dotted line is a fit obtained with a 1/L^ dependence. 
As it can be seen in this figure, the binding tends to zero in the infinite volume limit. This indicates 
that this two-fermion system has no bound state for this particular set of parameters. The situation is 
however the same for the whole range of parameters accessible in the numerical simulations. Although 
assuming that this pathology could be associated to the quenched approximation it is physically 
surprising that no any NN bound state could be generated unless the NN pairs creation (unquenched 
Yukawa model) are taken into account. 
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Figure 9: Two-body effective mass for a two- 
fermion 0+ state vs time with local source and 
fully smeared sources. L = 48, averaged over 
800 configurations, gi = 0.3, k = 0.118, and 
a/j, = 0.1. 




Figure 10: Binding energy versus lattice size L, 
for gL = 0.3, k = 0.118, and a/i = 0.1 aver- 
aged over 4000 samples for L = 16, • • • , 32, 2000 
for L = 36, and 800 for L = 48. Dotted line 
corresponds to a fit. 



4.3 Scattering lengths 

Since there is no bound state below the critical coupling constant, we can only access to the scattering 
properties of the state. The scattering observables cannot be obtained in Euclidean time in the infinite 
volume limit [28] but can be extracted from the volume dependent binding energy measured on finite 
lattices, like for instance the one plotted in figure [TOl The underlying formalism was developed by 
Luscher in [29\ I30j who gave a 1/L expansion of the the binding energy . In its leading order it reads: 

^ = Ip^^ (39) 

Taking the binding energy values of figure [TOl and equation (I39p . the NN scattering lengths ao have 
been extracted. The results corresponding to gi = 0.3, k = 0.118, and a/i = 0.1 are displayed in fig 
[TTl They manifest a constant behavior of oq a function of the lattice size L, indicating that the 
first order expression we used to extract the scattering length is justified. The dimensionless coupling 
constant of the nonrelativistc model for the parameters of figure [11] is G = 0.193. The corresponding 
non relativistic scattering length value, given by figure [5] is ^0 = ^o^ = —0.214 (solid line), quite 
close to its Born approximation (dashed line). 

This study has been performed for several values of gi- The dependence of ao on the coupling 
constant G is plotted in figure [121 foi' a lattice size of Lafi = 2.4 (L = 24, o/i = 0.1). One can see 
that the lattice results notably departs from the non relativisitc ones (solid line) and are above the 
Born approximation (dashed line). In figure [12] the NR scattering length has been computed using 
the lattice discretized potential of figure [3] As one can see, the regularization of the potential at the 
origin has no effect in the scattering length (indistinguishable from the continuum result in the figure) . 
This is due to the fact that the scattering length is a zero energy observable and therefore is not very 
sensitive to the details of the interaction. From the preceding analysis we conclude that the repulsive 
effect shown by the lattice data is not related to the lattice potential discretization at short distances. 
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Figure 11: Scattering length extracted from Eq. 
(|39p as a function of the lattice volume Lan, for 
gi = 0.3, K = 0.118 and a/i = 0.1 i.e. G = 
0.193, extracted from the data in fig. [TOl The 
solid line indicates the non relativistic results 
and the dotted the Born approximation (|32p . 




Figure 12: Scattering length vs G for a lattice 
volume La/i = 2.4, 2000 samples for each point 
up to G = 0.3 and 5000 thereafter. The solid 
line indicates the continuum non-relativistic re- 
sult and the dotted one the Born approximation 
(|32p . Blue circles represent the NR result ob- 
tained using the lattice potential in figure [3l 



The values of the accessible coupling constants extend beyond the Born regime but are still far 
from the pole behavior corresponding to the appearance of the first bound state displayed in figure 
[8l The difference between the lattice and NR results may indicate strong repulsive corrections. These 
kind of corrections are already manifested in the bound state problem when solving the same Yukawa 
model both in Light Front |31j and Bethe-Salpeter j32j ladder equations. 

5 Conclusion 

We have considered the quantum field theory solution of the simplest nuclear Yukawa model consisting 
on two identical nucleons interacting via a scalar meson exchange. The choice of the scalar coupling 
with respect to pseudo scalar one was taken in order to optimize the appearance of two-body bound 
states we were interested in. 

The problem has been solved using the standard lattice techniques, based on the path integral for- 
mulation of the theory on a discretized space-time. The meson field has been described by a discretized 
Klein-Gordon Lagrangian without self-interacting X(p^ term and the Dirac- Wilson discretization was 
chosen for the fermion. 

The resulting model is fully relativistic and was solved by neglecting only the NN loops generated 
from the meson field in the so called quenched approximation. This simplification is physically justified 
by the heaviness of the nucleon and is anyway implicit in all nuclear models. 

The numerical simulations were performed along the physical line fi/m ~ 0.6 where fi and m denote 
respectively the meson and nucleon renormalized masses. The solutions were found only for coupling 
constants below some critical value 5^3. Above this value the ubiquitous presence of fermion zero 
modes made the problem numerically unsolvable. The addition of a pseudo scalar coupling term does 
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not make the problem simpler. The present situation does not allow to judge whether this problem 
is related to the particular fermion discretization used. However the same problem was described in 
the past to affect naive fermions [2l]. This seems to indicate that the use of Wilson fermions is not 
responsible for the problem. 

The range of the accessible coupling constants is below the threshold value for producing two- 
body bound states, which in the non relativistic potential approximation turns to be g ~ 3.7 and in 
the Bethe-Salpeter one sensibly larger. In the accessible region of g the NN scattering length was 
calculated using the Luscher procedure. The values found were in agreement with the non relativistic 
models for low coupling constants but show strong repulsive effects when increasing g. 

We conclude that the quenched approximation of the Yukawa model (scalar coupling) is not able 
to produce two-nucleon bound states. Although this pathology is manifested in the quenched approx- 
imation it suggests that the non-relativistic results based on one-boson exchange potentials have no 
direct counterpart in the quantum field theory approach. 
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